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An  optimal  procedure  is  established  for  the  reconstruction  of  the  angular 
object  distribution  in  a given  field  of  view  (FOV).  The  object  is  coherently 
illuminated  and  located  in  the  far  tone  of  the  receiving  aperture.  The  proce- 
dure is  "uniformly"  optimal  in  the  sense  of  minimizing  the  statistical  rms 
difference  between  the  object  distribution,  modeled  as  a random  function  of  the 
angular  coordinates  and  its  reconstructed  image,  for  each  direction  belonging 
to  the  FOV.  The  observable  complex  amplitude  distribution  of  the  field  on  the  . 
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aperture  is  due  in  the  general  case  not  only  to  the  incident  field  scattered  by 
the  object  but  also  to  background  disturbance,  or  "angular  noise",  randomly 
distributed  inside  and  outside  the  FOV,  and  is  affected  by  "measurement  noise" 
that  is  random  errors  introduced  in  measuring  the  aperture  field.  The  recon- 
struction algorithm  consists  of  summing  a truncated  series  of  special  func- 
tions — prolate  spheroidal  for  tne  linear  case  and  their  generalization  for  two 
dimensional  apertures  — weighted  by  appropriate  coefficients.  These  coeffi- 
cients depend  upon  the  observed  aperture  field  and  upon  the  relative  power 
densities  associated  with  the  object  field  and  the  various  types  of  noise.  The 
series  is  truncated  to  a number  of  terms  ("effective  degrees  of  freedom"  of 
the  image)  determined  through  an  information  theoretical  method:  each  term 
of  the  series,  suitably  ordered,  provides  an  information  gain  less  than  the 
preceding  one,  and  the  information  gain  goes  rapidly  to  zero.  The  relationship 
between  information  transfer  and  mean -squared  error  for  each  term  in  the 
image  series  is  established.  Numerical  examples  are  discussed. 
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Optimum  Coherent  Imaging  of  a Limited  Field  of  View 
in  the  Presence  of  Angular  and  Aperture  Noise 


I.  INTRODUCTION 

In  recent  years  a substantial  amount  of  work  has  been  done  on  the  problem  or 
image  reconstruction  from  knowledge  of  the  field  distribution  on  a receiving  aper- 
ture. The  incoherent  case  has  been  discussed  by  Lr,  1 who  concluded  that  a com- 
plete restoration  of  a sky  temperature  distribution  is  possible  in  principle,  despite 

the  finite  ?ize  of  the  aperture.  The  problem  was  reconsidered  in  detail  bv  Buck 
2 

and  Gustincic  who,  exploiting  the  properties  of  spheroidal  functions,  showed  that 

only  a limited  number  of  fum  cions  was  useful  in  representing  the  reconstructed 

3 

temperature  distribution.  Rushforth  and  Harris  considered  the  effect  of  noise 
superimposed  on  both  the  object  and  the  image  distributions  and  took  into  proper 
account  the  advantages  given  by  the  a priori  knowledge  of  the  limited  extension  of  the 
object  in  its  plane.  Toraldo  di  Francia4 discussed  the  general  question  of  the 
degrees  of  freedom  of  an  image  and  pointed  out  the  difference  between  the  coherent 
(Received  lor  publication  28  April  1976) 

1.  Lo,  Y.T.  (1961)  On  the  theoretical  limitations  of  a radio  telescope  in  determin- 

ing the  sky  temperature  distribution.  J.  Appl.  Phys.  32:2052-2054. 

2.  Buck,  G.  J.  and  Oistincic,  Jacob  T.  (1967)  Resolution  limitations  of  a finite 

aperture,  IEEE  Transactions  on  Antennas  and  Prop.  A P-15,  No.  3:376-381. 

3.  Rushforth,  C.  K.  and  Harris,  R.  \V.  (1968)  Restoration,  resolution  and  noise.  J. 

Opt.  Soc.  Am.  58,  No.  4:539-545. 

4.  Toraldo  di  Francia.  G.  (1969)  Degrees  of  freedom  of  an  image,  J.  Opt.  Soc. 

Am.  59.  No.  7:799-804. 
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and  the  incoherent  case.  Bendinelli,  et  al^  established  expressions  for  the  recon- 
struction coefficients  — in  the  presence  of  measurement  noise— that  minimize  in  a 
statistical  sense  the  mean  square  difference  between  the  object  and  its  reconstruc- 
tion, integrated  over  the  extension  of  the  object.  All  those  studies  are  limited  to 
the  two-dimensional  case  of  linear  apertures. 

In  this  paper  the  question  of  image  distribution  reconstruction  *or  coherent 
illumination  is  reconsidered  from  first  principles,  without  any  a priori  assumption 
of  a particular  processing  system  behind  the  receiving  aperture,  that  is  without 
postulating  any  particular  "optical  spread  function. " The  object  is  assumed  to  be 
located  in  the  aperture  far  zone.  This  being  the  case,  the  set  of  spherical  waves 
constituting  the  scattering  contributions  from  each  point  of  the  obiert  can  be  con- 
sidered locally  planar  at  the  aperture.  The  image  reconstruction  is  therefore 
equivalent  to  the  determination  of  the  function  of  direction  characte-izing  the  com- 
plex amplitudes— referred  to  the  aperture  center— of  the  Plane  Wave  Spectrum 
(PWS)  into  which  the  aperture  field  due  to  the  object  can  t’  decomposed.  Because 
of  its  meaning  the  PWS  associated  with  the  object  will  be  denoted  as  object  angular 
distribution  or  simply  object  function,  and  mathematically  modeled  as  a complex 
random  function  of  the  angular  coordinates.  The  complex  coefficients  for  two  plane 
waves  incident  from  two  different  directions  are  assumed  to  be  statis'ically  uncor- 
related. However,  their  relative  phase  relationship  is  fixed  in  time.  This  is 
equivalent  to  considering  the  coherent  case  only.  It  is  assumed  that  the  object  dis- 
tribution of  interest  is  angularly  limited  to  a certain  a priori  assigned  field  of  view 
(FOV1.  We  want  to  reconstruct  its  complex  values  by  linear  processing  of  the  com- 
plex  amplitude  of  the  field  distribution  observed  at  the  receiving  aperture.  The 
observed  values  are  due  not  only  to  the  incident  field  scattered  from  the  object  but 
also  to  background  interference,  or  "angular  noire",  generated  by  scatterers  from 
inside  and  outside  the  FOV  in  the  aperture  far  zone  and  by  "measurement  noise" 
locally  introduced  at  the  aperture  in  the  measurement  process. 

In  this  paper  a linear  reconstruction  procedure  will  be  established  by  using  the 
methods  of  statistical  estimation.  The  procedure  is  "uniformly  optimal"  in  the 
FOV  in  the  sense  of  minimizing  the  statistical  rms  difference  between  the  object 
distribution  and  its  reconstructed  image,  for  each  direction  of  interest  (ratner  than 

5.  Bendinelli,  M. , Consortini.  A.,  Ronchi,  L.,  and  Frieder.,  B.  R.  (19741  Degrees 
of  freedom  and  eigenfunctions  for  the  noisy  image,  J.  Opt.  Soc.  64,  No.  11: 
1498-1502.  

3 

’The  situation  is  similar  to  that  examined  by  Rushforth  and  Harris  for  a somewhat 
different  situation,  since  they  considered  the  optical  property  of  the  system  speci- 
fied a priori  through  a point  spread  function.  However,  unlike  the  case  considered 
in  Ref.  3 in  this  paper  all  the  scatterers  are  assumed  to  be  in  the  far  zone  of  the 
aperture.  Consequently,  the  background  noise  in  the  object  plane  in  Ref.  3 is 
replaced  by  a PWS  modeled  as  a random  function  of  direction  and  denoted  as 
"angular  noise. " 
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the  rms  reconstruction  error  integrated  over  the  FOV,  as  for  example  in  Ref.  51. 

The  prolate  spheroidal  functions,  for  the  linear  aperture,  and  their  generalizations 
for  two-dimensional  apertures,  play  a fundamental  role  in  the  analysis,  as  they  did 
in  most  of  the  other  related  work.  In  previous  work,  however,  the  use  of  spheroidal 
functions  stemmed  from  their  being  eigunfunctions  of  the  integral  equation -defining 
the  imaging  operation-whose  kernel  was  the  optical  spread  function.  Their  rele- 
vance to  reconstruction  is,  in  a sense,  far  more  fundamental.  This  was  recognized 
4 2 

byToraldo  di  Francia  and  Buck  and  Gustincic.  The  appearance  of  spheroidal 
functions  and  their  two-dimensional  generalizations  (discussed  in  this  paper  for 
relatively  arbitrary  geometries),  is  shown  to  be  a natural  consequence  of  the  struc- 
ture of  the  inhomogeneous  integral  equation  defining  the  statistically  optimal  linear 
processing  of  the  observed  aperture  field.  The  reconstruction  algorithm  consists  of 
summing  a truncated  series  of  functions  (prolate  spheroidal  or  their  generalization) 
weighted  by  coefficients  depending  upon  the  observed  aperture  field  and  upon  the 
statistical  second  moments  of  the  object  distribution  process  and  of  the  various  types 
of  noise.  The  series  is  truncated  to  a number  M of  terms  ("Effective  Degrees  of 
Freedom")  determined  through  considerations  of  information  theory.  Each  term  of 
the  series,  suitably  ordered,  provides  an  information  gain  less  than  the  preceding 
one.  The  number  M is  such  that  no  advantage  is  obtained  by  adding  additional  terms 
to  the  serfes  for  the  image  reconstruction.  For  each  additional  term  the  uncondi- 
tional entropy  and  the  entropy  conditioned  to  the  presence  of  a given  random  scene 
become  asymptotically  equal.  Hence  the  information  gain  tends  to  zero.  In  terms 
of  rms  errors  integrated  over  the  FOV  this  fact  means  the  following:  the  integrated 
rms  error  associated  with  a term  of  order  i > M in  the  reconstruction  series  is 
essentially  equal  to  the  variance  of  the  corresponding  terms  of  the  expansion  of  th* 
random  scene  integrated  over  the  FOV.  Therefore,  those  terms  in  the  series  con- 
tribute only  to  noise  in  the  reconstructed  image.  The  number  M is  proportional  to 
the  FGV-apenure  product,  unless  very  low  and  possibly  unrealistic  values  of  dis- 
turbance are  present.  In  such  cases  the  possibility  exists  of  an  improvement  of  the 
reconstruction  accuracy  beyond  the  limits  suggested  by  the  classical  optical  theory. 
The  domain  of  application  of  the  reconstruction  method  here  proposed  is  re- 
stricted to  microwave  frequencies.  In  fact,  the  assumption  that  the  field  can  be 
observed  at  all  points  of  the  aperture  implies  that  the  field  is  a classical  electro- 
magnetic field,  which  requires  a large  number  of  quanta  per  degree  of  freedom. 

This  condition  is  encountered  at  optical  frequencies  only  with  extremely  intense 
fields. 
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2. 


MATHEMATICAL  MODEL  AND  PRELIMINARY  RESULTS 


2. 1 Object  Function  and  Angular  Noiw 

We  will  consider  apertures  and  FOV's  which  are  generally  two-dimensional. 

It  will  be  straightforward  to  simplify  notation  and  results  in  order  to  deal  with  the 
simpler  case  of  linear  aperture. 

Let  A be  a receiving  aperture  on  the  x,y  plane  (whose  area  will  also  be  indi- 
cated by  the  same  letter  A).  It  is  expedient  for  mathematical  reasons  to  assume 
that  the  aperture  has  a point  symmetry  with  respect  to  its  center.  This  mear  s,  if 
the  origin  is  assumed  coincident  with  the  aperture  center  and  an  aperture  point  is 
located  at  <x,y),  that  there  exists  another  aperture  point  at  (-x, -y).  Except  for 
this  constraint  of  point  symmetry,  the  geometry  is  arbitrary. 

Let  a position  vector  on  the  aperture  plane  be  given  by: 

x = x x + y y , (1) 


x and  y being  unit  vectors  in  x and  y directions.  The  rectangular  coordinates  in 
the  wavenumber  plane  are  related  to  the  angular  coordinates  of  a standard  spherical 
system  as  follows: 


2ir  . a 

u = -j-  sin  8 cos  tp. 

(2) 

2*  . . . 

v = — sin  8 sin  , 

(3) 

where  X is  the  wavelength.  A position  vector  ir.  the  u, 
duced: 

v plane  is  conveniently  intro- 

u = u x -f  v Y . 

(4) 

Let  A(x)  be  a function  equal  to  unity  in  the  points  of  the  aperture  and  zero  elsewhere. 
/ I forx<A 

A(x)  = J , (5) 

* 0 for  *4  A 

the  ’’aperture  function"  is  defined  as  follows 

♦ 4-  « 

«<u)  = / J A(x)ej-id2x  . 

-«e  -«c 

(6) 
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where,  of  course. 


u-  x = u x + v Y , 

2 

and  d x is  the  element  of  the  area  in  the  x,y  plane.  The  field  of  view  is  defined  as 
the  domain  S of  the  u,  v plane,  inside  which  is  located  the  object  distribution  to 
^e  reconstructed.  S is  assumed  to  have,  like  A,  a point  symmetry  about  its  center 
in  the  u,  v plane.  The  same  letter  S will  denote  the  area  of  the  FOV  in  the  wave- 
number  plane.  Also  it  proves  convenient  to  introduce  the  function: 

( 1 for  ufS 

S(u)  = ~ (7) 

I 0 for  u 4 S 

Let  the  object  angular  distribution  be  a random  function  g (u),  different  from 
zero  only  in  S,  which  will  be  called  the  object  function.  The  function  gQ(ui  is  the 
complex  amplitude  (referred  to  the  aperture  center)  of  the  PWS  representing  the 
object.  As  mentioned  in  Section  1,  we  assume  that  the  values  of  g (u)  for  two  differ- 
ent arguments  are  statistically  uncorrelated.  This  physically  means,  intuitively 
speaking,  that  in  the  FOV  we  have  no  a priori  information  of  how  the  object  function 
in  a certain  direction  u of  S affects  probabilistically  the  value  of  the  object  function 
in  a neighboring  direction.  Thus,  if  we  denote  the  statistical  average  operator  by 
E. 


K [ S(u)  go<u>  S(v)  go<v>  j = 4*2  «<u  - v>  P S(u>  , (8) 

where  v,  like  u,  represents  an  arbitrary  direction  and  6(u>  is  the  two-dimensional 
— — 2 — 
impulse  function.  Pj.  has  the  physical  meaning  of  4*  times  the  power  per  unit 

aperture  area  incident  from  an  unit  area  of  wavenumber  plane. 

A random  function  Nj(u>,  different  from  zero  only  in  the  domain  S of  the  u,  v 

plane,  is  introduced  in  order  to  represent  the  PWS  associated  with  background 

coherent  disturbance,  or  angular  noise,  incident  from  directions  belonging  to  the 

FOV.  Nj<u).  like  gQ(u),  is  an  uncorrelated  random  function  homogeneous  in  S. 

Hence: 


E ^S(u)  Nj(u)  S(v)  N j (v)  j = 4*2  P,  S(u>  6(u  - v>  . (9) 

Noise  is  present  outside  the  FOV.  It  may  simply  mean  that  the  part  of  the  scene 
outside  the  FOV  is  of  no  interest  and  is  therefore  considered  a disturbance.  Again 
such  a noise  Njfu)  is  uncorrelated: 


[ 1 - S(u)j  N2(u)  1 - S(v)  j 


N2(v>  | 


= 4w2P. 


[l  - S(u)j6(u  - v)  . (10) 


and  a reasonable  assumption  is  that  the  statistical  cross  correlation  between  go<u), 
Nj(u)  is  zero  for  all  values  of  the  arguments. 

In  the  analytical  model  expressed  by  (10)  the  angular  noise  extends  throughout 
the  whole  wavenumber  plane,  although  of  course  Ngfu)  should  be  zero  outside  the 
circle  of  the  wavenumber  plane  corresponding  to  real  directions,  tnat  is,  outside 
the  "vis.ble  space"  G,  defined  as  the  set  of  points  of  the  wavenumbers  plane  such 
that 


1 ul  < 2W > . 


The  assumption  (10)  however  is  useful  in  simplifying  the  subsequent  development, 
and  has  been  adopted  in  most  of  the  previous  work  in  this  area.  In  Appendix  C,  it 
will  be  shown  that  negligible  error  is  committed  by  assuming  the  validity  of  (10), 
because  in  the  wavenumber  plane  the  visible  space  has  a much  greater  extension 
than  the  FOV  we  are  looking  at. 


2.2  Aperture  Field  and  Aperture  Noi*e 

The  received  aperture  field  is  the  superposition  of  the  plane  wave  spectra 
representing  the  object  and  the  angular  noise.  Therefore  in  our  scalar  approxima- 
tion it  takes  the  form  of  a Fourier  Transform: 


•* « 


[■ 


N2<u> 


j e'iH'id2. 


(ID 


which  defines  f(x)  as  a random  function.  If  ( , like  x represents  a point  of  the  aper- 
ture from  (9  to  11),  it  easily  follows  that: 

+»  ♦ » . 

E [f(x*  f <{>J  * A(x)  i.i|)  j J Js<u>  (Ps  + Pj ) + [l  - S(u)j  P2  J ej-(x  ' * *d2u, 

(12) 


that  is. 


K [f(x)  f (§)]  = A(x)  A({>  ^ -ly  Os(x-l)  ( PS  + Pi  ' P2)  + P2  8<-  'i’j  • 


(13) 


where  the  characteristic  function  of  the  FOV  has  been  introduced: 


10 


+ 00 


(14) 


i 


f 

l 


+ 00 

(TS(*>  = f J S(u)e  — d^u  . 

-00  -00 

To  complete  our  model  let  us  assume  that  additive  noise  n(x)  is  present  on  the 
aperture.  This  means  that  "measurement  errors"  are  present  in  the  observation 
of  the  aperture  field  distribution.  In  practical  cases  this  noise  may  be  of  a multi- 
plicative nature.  However,  the  additive  model  simplifies  the  mathematics  and  is 
close  to  reality  in  many  cases.  Thus  .ne  "observable"  field  dist-ibution  on  the 
aperture  is 

o(x)  = fix)  4 n(x)  , (15) 

a random  function  having  zero  mean,  and  we  assume-  that  the  measurement  noise  is 
uncorrelated: 

E [ A(x)  A(£)  nix)  n*(|>]  = rr * A(x)  6(x  -^>  . (16) 

This  assumption  is  made  customarily  for  mathematical  convenience,  5 and  leads 
to  the  well  known  difficulty  of  implying  an  infinite  noise  power  for  a continuous 
aperture.  The  difficulty  is  merely  formal  however  and  disappears  if  the  rms  com- 
ponents of  aperture  noise  with  respect  to  a system  of  orthogonal  functions  are  con- 
sidered (see  Appendix  C).  Since  the  object  function  and  the  noise  are  uncorrelated 
one  easily  obtains  by  recoiling  (9  to  11)  and  (15): 


E ^o(x)  o (|)J 

I 

= A(x)A({) 


4» 


'-7  «s(i  -i>  (PS  4 P1  * P2  ) 4 (P2  + °n  ) »<*  -i>  • < 


17) 


3.  LINEAR  ESTIMATION  OF  THE  OBJECT  FUNCTION 
3.1  Sutemrat  of  the  Problem 

We  want  to  estimate  the  object  function  gQ(u)  by  linear  processing  of  the  observ- 
able aperture  field  distribution  o(x).  Hence  the  estimator  $(u)  will  have  the  general 
form  of  linear  transformation: 

$(u)  =,  S(u)  f f H(u,  x)  o(x)  d^x  , (18) 
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where  o(x)  is  given  by  (15)  and  the  kernel  H(u,x)  is  to  be  determined.  To  accom- 
plish such  a task  we  require  that  for  any  direction  u belonging  to  the  FOV,  the 
statistical  average  of  the  rms  estimation  error  be  minimum: 

2 

e2  = E | |(u)  - gQ(u)  | = min.  (19 

As  shown  in  Appendix  A this  leads  to  the  integral  equation: 

S(u)  A(x)  ej-  - 


/ P.  - P9\  i , , 9 + 

:u)A(x)  (l + _ip— il-ij-  //  os(x  - J)  H(u,j)  d £ + H(u,  x)  . 

' S'  A S 


= S(u) 


whose  solution  for  H(u,  x)  will  provide  the  optimum  linear  processing  scheme. 


3.2  Optimum  ftocemor 

To  solve  Eq.  (20)  consider  the  associate  homogeneous  eigenvalue  equation: 

\VX)  ‘7?  //  * Vl)d2i-  (21 

4ir  A 

The  kernel  ts  symmetric  real  and  positive  definite,  therefore  the  eigenvalues  are 
positive  numbers.  The  eigenfunctions  can  be  chosen  to  be  real  and  constitute  an 
infinite  dimensional  set  complete  and  orthogonal  in  A.  They  will  be  ordered  for 
increasing  values  of  X.,  starting  from  i = 0.  We  will  normalize  them  to  unity: 

;/a  *,(«>  v«>  d2. . «lk . <22 


H(u,  x)  V a <u)  <Mx) 
cln  * - * — 


and  insert  this  expression  into  (20),  obtaining 
S(u)  A(x)  ej- - 


I / P,  - P«\  * P„  ♦ <r2  » l 

= S(u)  A(x)  ll  Z1  V a (u)  A.  <Mx>  ♦ p P £ Vu) 

~ \ *S  ’ i= 0 1 ” 1 1 ~ *S  i 0 1 “ 1 ” J 
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To  find  the  coefficients  a.(u)  we  need  to  introduce  the  set  of  functions: 


*.(u)  = //  <Mx)  e-*—  -d2x  . 
i - J JA  i - 


(25) 


which  because  of  Parseval's  theorem  form  an  orthogonal  (although  not  complete) 
set  on  the  wavenumber  plane: 


f " ,+0°  *o  2 

/ / ^(u)  d £ = 4ir  6.^ 


(26) 


-00  -00 


A brief  discussion  of  the  properties  of  the  set  {♦.(u))  is  given  in  Appendix  B.  Now 
multiply  both  the  sides  oi  (24)  by  (ik(x)  and  integrate  the  aperture  A.  By  using 
orthogonality  and  the  definition  (25): 

r p _ p P + 

S(u)  *k(u)  = S(u)  (l  + — !p 1)  ak(u)  Ak  + — - ak(u)  (27) 


and  thus 


ak(u)  = 


S(u)  * (u) 

r-V)v 


P_  + a2 

2 n 


(28) 


The  function  H(u,  x)  therefore  is  from  (23): 


H(u,  x) 


(29) 


Hence  the  estimator  of  the  object  function  is  found  from  (18)  and  (29)  to  be: 


g(u)  = S(u)  ^ 


*k(u> 


k =0/  P,  - P2\ 


P2+6-  A 


V jf  V-  0<-)  d2x  • 


(30) 


completing  the  formal  solution.  The  estimate  (30)  will  be  called  the  best  image  or 
simply  the  image  of  the  object  function  gQ(u). 

In  any  practical  reconstruction  algorithm  the  series  (30)  is  truncated  after  a 
limited  number  M of  terms.  In  the  next  section  the  determination  of  the  number  M 
will  be  treated  by  considering  the  information  transfer  from  object  to  image 
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associated  with  each  term  of  (30).  It  will  be  shown  that  unless  measurement  and/or 
angular  noise  have  extremely  low  values,  essentially  no  advantage  is  obtained  by 
considering  terms  of  order  greater  than  N where  N is  defined  as 

2 

N = smallest  integers  AS/4ir  . (31) 

Such  an  information  theoretical  approach  has  been  preferred  to  the  alternative  ap- 
proach consisting  of  minimizing  in  some  sense  the  reconstruction  error.  In  fact 
the  1.  Ur  method  leads  to  certain  mathematical  convergence  difficulties  as  a 
consequence  of  the  assumption  of  a flat  spectral  density  for  gQ(u),  as  discussed 
briefly  in  Section  5.  Also  the  method  used  here  provides  an  insight  into  the 
mechanism  of  information  transfer  from  object  to  image. 


4.  INFORMATION  TRANSFER  FROM  OBJECT  TO  IMAGE 


4.1  Ini'ormalion  Gain  Aiiociated  with  the  Observable  Quantities 

In  the  processing  scheme  described  in  the  previous  section,  the  quantities 


O.  = j'f  *k(x)  o(x)  d2x 
A 


(32) 


are  the  "obseivables, " and  are  random  variables  whose  statistical  properties  are 
induced  by  these  of  the  random  function  o(x).  Since  o(x)  has  zero  mean,  we  have 
also 


E[0.]  = 0 . 


The  second  order  moments  of  the  observables  (32)  are  found  by  using  the  ortho- 
gonal properties  of  the  set  of  functions  {^(x))  , and  the  expression  (17).  In  this 
way  one  obtains: 


O.  O. 

t k 


x p ♦ x. 

i rS  i 


V 


(1 


P2+*n 


6, 


ik 


(33) 


Thus  the  ob«»rvables  are  uncorrelated.  Also  because  of  their  definition  (32),  they 
are  integrals  of  random  functions.  Thus  by  invoking  the  central  limit  theorem,  we 
may  assume  that  their  probability  densities  are  normal.  The  same  argument  can 
be  made  in  regard  to  the  conditional  probability  density  assuming  a certain  object 
function  gQ(u).  For  the  unconditional  probability  density  the  variance  is  from  (33): 


14 


IT 


(T2  = X.  Pc  + X.  P,  + (1  - X.)  P,  + 62  . 
l i b i 1 i i n 


and,  for  the  conditional  probability  density  assuming  gQ(u),  the  variance  is: 

cr2  = X.  P.  + (1  - X.)  P04  ct2  . (35) 

lg  i 1 l i n 

that  is  the  variance  of  O.,  due  only  to  the  various  types  of  noise. 

For  each  observable  O.,  the  entropies  associated  with  the  unconditional  and 

conditional  probability  densities  are  found  under  the  assumption  of  normal  distribu- 

6 

tions.  Thus,  for  the  unconditional  entropy  we  have,  by  using  here  base  2 log- 
arithms: 


H(i)  = log 2 (<ri  yfT*)  * 2~T5g~2~ 


with  <T.  given  by  (33).  For  the  conditional  entropy  assuming  gQ(u)  we  have  similarly: 


hbih  = ioe2  fasl  V7.1 


with  rtg.  given  by  (35).  It  is  known  that  for  a continuous  probability  density  the 
entropy  can  be  defined  to  the  extent  of  an  arbitrary  additive  constant.  This,  how- 
ever, does  not  create  any  problem  if  one  deals  with  a difference  of  two  entropies. 
The  information  gain  associated  with  the  observable  O.,  is  in  fact  equal  to  the  dif- 
ference of  the  two  entropies  (36)  and  (37),  that  is. 


I = H(i>  — H (i)  = logg  ~- 

gi 


or.  recalling  (34-35): 


Ij  = log2 


P2  + "n 


Pi'F2  . . fV"n" 


6.  Woodward,  P.  M.  (1953)  Probability  and  Information  Theory,  with  Application 
to  Radar  Per  gam  on  Press,  London,  chap.  1,  pp  21  -2b. 
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Since  the  A.'s  are  a monotonic  sequence  decreasing  with  i,  with  a step  behavior 
when  i reaches  the  value  N given  by  the  aperture  - FOV  products  (31)  (Appendix  B) 
we  have: 

lim  I.  = 0 , (40) 

i -*  oc 


and  an  asymptotic  expression  for  I.  is  found  to  be: 
i Pc  A- 

1:  -LT  log2e-  (41) 

' 2 


Equation  (41)  shows  that  the  information  gain  is  asymptotically  proportional  to  the 
characteristic  number  A.  associated  with  the  ith  degree  of  freedom  multiplied  by  a 
suitably  defined  signal  to  noise  ratio.  Practically,  unless  the  signal  to  noise  ratio 
is  very  large,  the  effective  number  of  degrees  of  freedom  of  the  image,  that  is, 
the  number  M of  the  reconstruction  terms,  can  be  chosen  equal  to  N in  (31),  be- 
cause terms  of  higher  order  carry  essentially  no  information.  This  conclusion  is 
corroborated  by  the  reconstruction  error  analysis  in  tile  next  section,  and  by  numer 
ical  results  given  in  the  sequel. 


4.2  RMS  Recoiutniction  Error  Averaged  in  the  FOV. 

Informational  Compared  to  Statistical  Approach 

It  is  shown  in  Appendix  B that  the  set  of  functions  {♦.(u)}  is  orthogonal  and 
complete  in  S.  Thus,  through  a standard  procedure,  we  can  represent  the  object 
function  by  the  expansion: 


g0<u>  = N_— J 
i=D  4ir  A. 


*i'i»  J/„ 


i d2u  • 


(42) 


Equation  (42)  is  obtained  without  resorting  to  any  statistical  consideration.  However 
it  could  have  been  obtained,  less  directly,  through  an  application  of  Parseval’s 
theorem  to  the  integrals  in  Eq.  (30)  for  the  estimate,  in  the  limiting  case  of  absence 
of  noise: 


This  means  that  in  the  ideal  conditions  of  Eq.  (43)  the  object  function  can  be  exactly 
reconstructed  in  principle,  despite  the  limited  size  of  the  aperture.  This  result  is 
perhaps  surprising,  but  well  established,  and  a brief  discussion  of  this  point  is  given 
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in  Ref.  4.  It  is  not  of  practical  interest  because,  as  soon  as  any  random  disturb- 
ance intervenes,  a very  different  conclusion  is  reached.  In  this  section  we  will 
consider  the  rms  errors  associated  with  each  term  in  the  estimate  (30),  and  we  will 
show  that  for  terms  for  which  the  observables  (32)  carry  zero  information,  the 
reconstruction  error  is  equal  to  the  a priori  variance  of  the  randor.  object  distribu- 
tion. The  result  shows  the  equivalence  of  the  informational  and  statistical  approach 
in  establishing  the  number  of  terms  necessary  in  the  reconstruction  procedure  (30). 

Consider  the  rms  reconstruction  error  averaged  in  the  FOV,  defined 


e2=  E| 


If  80M  - g(u) 


d2u 


(44) 


In  order  to  avoid  complicated  and  after  all  unnecessary  convergence  difficulties,  we 
assume  for  the  moment  that  gQ(u)  is  represented  satisfactory  by  the  series  (42) 
truncated  to  a finite  number  of  terms.  By  using  (B4)  of  Appendix  B,  (30)  and  (42), 
we  can  write  (44)  as  : 


l 4tr2  A.  E 

//a  Vj(*>  o(x)  d2x 

liS*i(“>  «o(H>d2*i 

(l  * ?1  " P2)\  , ?2  "" 

4ff2x 

l 

V ps  / i ' ps 

2 

1.  (45) 


where  the  number  of  terms  of  the  series  for  the  time  being  is  left  unspecified.  By 
applying  Parseval's  theorem  to  the  first  of  the  integrals  in  each  term  of  the  series 
(45)  and  recalling  that  the  random  object  function  is  uncorrelated  with  the  disturb- 
ance, one  obtains  through  simple  manipulations: 


W . P2 * °n 

B + p 

e 2 = 4e2P,V_^ l IL- 

a S 


i / P,  - P-\  P,  4 a* 
(1  + — P T~i  V "‘pc 


(46) 


By  recalling  (39),  the  expression  of  the  rms  error  (46)  turns  out  to  be 


(47) 


According  to  the  discussion  of  Section  4.  1,  the  information  associated  with  the  higher 
•«*erms  of  the  reconstruction  series  is  essentially  zero.  For  such  terms  Eq.  (47) 
consistently  shows  that  the  reconstruction  errors  tend  to  a constant,  proportional 
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to  the  density  of  power  per  aperture  unit  area  and  per  unit  area  of  wavenumber 
plane  associated  with  the  object  function.  A further  insight  on  the  meaning  of  Eq. 
(47)  is  obtained  by  noting  that  the  rms  of  any  term  of  the  series  Eq.  (42)  for  the 
object  function  integrated  in  the  FOV  is: 


//  J3SJ4  A 

//  «i<v)  g (v)  d2v 

2 " 

S (4ir  X.) 

S 

(48) 


independent  of  the  index  i.  It  can  be  established  that  Eq.  (48)  is  also  equal  to  the 

asymptotic  value  of  each  term  of  Eq.  (47)  when  the  index  i tends  to  infinity  (and  I 

therefore  tends  to  zero).  This  means  that  the  integrated  rms  errors  associated 

with  higher  order  terms  in  the  reconstruction  series  are  equal  to  the  a priori 

integrated  mean  squared  value  of  the  object  function  itself.  Thus  adding  terms 

carrying  zero  information  does  not  improve  or  deteriorate  the  average  mean 

squared  error  of  the  reconstructed  image.  Again  we  reach  the  conclusion  that,  not 

unexpectedly,  the  information  theoretical  and  statistical  appro,  ches  are  equivalent 

in  expressing  the  circumstance  that  the  uncertainty  associated  with  higher  order 

> 

reconstruction  coefficients  is  equal  to  the  a priori  uncertainty  of  the  corresponding 
terms  in  the  representation  of  the  random  scene.  These  terms  are  therefore 
useless  in  the  reconstruction  procedure. 


5.  PARTICULAR  CASES 

Two  important  particular  situations  will  be  now  considered: 

(a)  Imaging  a Limited  Sector  of  a Wide-Angle  Random  Scene 
In  this  case  we  assume  that  the  aperture  receives  energy  from  a statistically 
homogeneous  random  scene,  only  a limited  part  of  which  - the  one  belonging  to  the 
FOV  - is  actually  of  interest.  This  case  can  be  modeled  by  assuming  that  gQ(u)  and 
Njlu)  are  different  parts  - inside  and  outside  the  FOV  - of  the  same  random  process, 
so  to  speak.  Consequently  we  assume  in  our  model  that  Nj(u)  is  identically  zero. 
Thus: 


Pl-  ° : PS=  P2‘ 


The  reconstruction  procedure,  Eq.  (30),  is  in  this  case: 

,2 

‘ £ Vjt'U'  J ) 'JW  r 

'n 


|(u) 


u)  = S(u) L-j-  T Vk'u>  //  *.(x)  O(x)  d2x  , 

ff  k A 


1 ♦ 


(49) 


(50) 
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Equation  (50)  shows  that  the  reconstruction  algorithm  docs  not  depend  upon  the 
2 

ratio  Pc/fT  , which  appears  in  a scale  factor.  In  other  words  the  form  of  the 
o n 

estimator  g(u)  is  independent  of  the  severity  of  the  measurement  errors  and  of  the 
power  received  per  square  meter  af  aperture.  However  the  measurement  error 
affects  the  information  transfer  from  the  scene  to  the  image,  which  for  the  ith 
reconstruction  term  is  from  Eq.  (35): 


I. 

t 


(51) 


In  the  next  section  we  will  see  that  for  this  case  the  optimal  reconstruction  proces- 
sor is  equivalent  to  an  isoplanatic  optical  system. 

(b)  Objects  in  the  FOV  Immersed  in  a Wide-Angle  Background  of 
Angular  Noise 

Assume: 


(52) 


The  reconstruction  algorithm  is: 


v *’i(£)  rr  2 

g(u)  = S(u)  ' «-  If  *.(x)  O(x)  d x , 

~ t P„  ♦ n_  A 1 “ ~ ~ 


r 

V” F 


(53) 


which  does  depend  upon  the  level  of  background  noise  and  measurement  noise  with 
respect  to  the  object  power.  The  information  transfer  is  from  (39): 


I.  = 

i 


X + 

l 


i loe2 


pc  + V 


(54) 


and  goes  rapidly  to  zero  for  i > N. 

A numerical  example  is  provided  to  clarify  the  meaning  of  the  theoretical  re- 
sults obtained. 

Example: 

Consider  the  one-dimensional  case  of  a linear  aperture  of  length  2a.  The  angu- 
lar extent  of  the  FOV  is  26  , thus  in  the  wavenumber  axis  its  extent  is: 

s 
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“ “ - — — r — jiu  C? 

s X s 


The  one  dimensional  equivalent  of  (31)  is: 

2a  2 u 


N = smallest  integer  > 


TT 


2 sin  0 
X s 


Put: 


and  assume  for  our  numerical  example: 
c = 8,  that  is,  N = 5 , 

The  characteristic  Eq,  (21)  in  the  monodimensional  case  becomes  the  eigenvalue 
equation  for  prolate  spheroidal  functions,  with  parameter  c.  The  first  eight  eigen- 
values for  i = 8 are  listed  in  Table  1. 

The  information  transfers  associated  with  the  various  terms  in  the  reconstruc- 
tion series  for  cases  (a)  and  (b)  considered  above  are  now  numerically  evaluated. 
Case  (a). 

2 

For  various  values  of  the  ratio  Pc/o  the  information  contents  i.  of  the  various 

on  i 

terms  in  the  image  reconstruction  are  listed  in  Table  2.  if  is  apparent  that  the  F's 
for  i > 4 decrease  very  rapidly.  In  fact  I,.  is  already  negligible.  Notice  the  weak 
dependence  of  the  higher  terms  of  the  image  series  upon  the  signal  to  noise  ratio. 

As  expected  from  the  theoretical  considerations,  five  terms  only  carry  practically 
the  total  information  available. 

Case  (bl 

In  Table  3 the  reconstruction  coefficients  for  case  (b)  have  been  Dated,  normal- 
ized to  the  first  (i  = 0)  for  various  values  of  the  appropriately  defined  signal  to 
noise  ratio.  The  information  transfers  are  listed  in  Table  4.  Again  we  reach  the 
conclusion  that  N = 5 terms  are  sufficient  for  the  optimum  reconstruction  for  small 
or  moderate  signal  to  noise  ratios.  However  for  a signal  to  noise  ratio  of  the  order 
of  30  d3  or  greater  the  possibility  of  a moderate  increase  in  the  accuracy  of  image 
reconstruction  emerges.  An  inspection  of  Table  4 shows  that  Ig,  Ig,  and  i^  are  not 
negligible  when  the  signal  to  noise  ratio  is  equal  to  10^. 
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6.  DISCUSSION  AND  FINAL  REMARKS 


As  customary,  in  our  treatment  the  scalar  approximation  has  been  assumed. 
For  limited  FOV  - ten  degrees  or  so  - it  is  believed  that  no  substantial  error  is 
committed.  On  the  other  hand  the  mathematical  complexity  is  substantially  re- 
duced, because  the  analytical  machinery  of  double  orthogonal  functions  can  be 
applied. 

The  physical  reason  why  higher  terms  in  the  series  (30)  are  of  little  or  no  use 
in  the  reconstruction  algorithm  stems  from  the  fact  that  the  "energy"  associated 
with  the  functions  <l>.(u)  is  essentially  confined  within  the  FOV  for  i s N,  and  outside  for 
i > N (see  Appendix  B).  Thus  higher  order  reconstruction  terms  contribute  little 
to  a faithful  reproduction  of  the  random  scene  inside  the  FOV.  On  the  other  hand 
by  applying  Parseval's  theorem  to  the  integrals  in  (30)  it  is  seen  that  higher  order 
terms  are  those  maximally  affected  by  the  angular  noise  incident  from  outside  the 
FOV.  To  get  a better  heuristic  insight  of  the  nature  of  the  reconstruction  algorithm 
consider  the  case  (a)  of  the  preceding  section.  In  this  case  we  singled  out  the 
sector  S - the  FOV  - from  a wide  angle  random  scene.  Truncating  the  appropriate 
estimator  (50)  to  N terms,  (N  being  the  smallest  integer  greater  th  .r.  the  aperture  - 
FOV  product),  or  letting  the  sum  be  extended  to  infinity  yields  informationally  equi- 
valent expressions  as  follows  from  the  discussion  of  Section  4.  We  invoke  now  the 
representation  for  the  aperture  function  (easy  to  establish): 

00 

..•j-v)  = -lli  *.<u)**(v)  (55) 

' 4ir  i=0  1 “ 1 “ 

valid  for  any  u,  v,  and  uniformly  convergent  for  u and  v both  belonging  to  S.  7 By 
taking  the  upper  limit  of  the  sum  in  (50)  equal  to  infinity,  by  applying  Parseval’s 
theorem  to  the  various  integrals,  and  interchanging  the  sum  with  the  integral,  one 
gets  (neglecting  a constant  factor): 

+ 0C  4 *° 

|(u)  - S(u)  — / J »A<u  ' d2Jl  • <r,6) 

4t  -oc 


where 

O(u)  = [gQ(u)  ♦ N j (u)  j S(u)  4 N,(u)  [l  - Sfulj 

4 //  n(x)  e-*  — ’ — d2x  (57) 

A 

7.  Hildebrand,  Francis  B.  (1965)  Methods  of  Applied  Mathematics,  2nd  Edition, 
Prentice  Hall,  N.J.,  pp  314-TT5: 
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is  the  angular  distribution  of  the  incident  field  plus  an  equivalent  angular  measure- 
ment noise— the  last  term  in  (57).  Equation  (56)  has  a simple  physical  interpreta- 

2 

tion.  For  case  (a)  the  optimal  processing  is  independent  of  the  ratio  Pj,/on  and  is 
equivalent  to  the  operation  performed  by  a clear,  aberration  free  (isoplanatic) 
optical  system,  with  "point  spread  function"  equal  to  the  aperture  function  rr^(u), 
the  observed  output  being  limited  to  directions  belonging  to  the  FOV. 

In  different  cases,  like  (b)  in  the  previous  section,  the  optimal  processor  de- 
pends upon  the  a priori  knowledge  of  signal-to-noise  ratio.  However,  because  of 
the  cutoff  properties  of  A.,  it  can  be  inferred  that,  only  when  the  noise  is  very 
small,  an  improved  restoration  is  obtained  by  using  the  optimum  processor  instead 
of  the  optical  operation  (56). 

To  conclude  these  remarks  we  want  to  establish  an  approximate  but  very  inter- 
esting expression  for  the  total  information  that  can  be  extracted  from  a limited  FOV 
in  the  presence  of  angular  and  measurement  noise.  By  invoking  the  cutoff  property 
of  the  characteristic  number  A.,  we  assume  that  the  actual  values  A.  can  be  approxi- 
mately replaced  by  1 for  i s N,  and  by  0 for  i > N.  Also  assume  [see  (31)]: 


Then  from  (39)  we  obtain  fcr  the  total  information  1: 


. ~ AS 

* 7"? 

4 w 


i log2  (* + 


P9  + ffn 

l n 


which  closely  resembles  the  classical  expression  for  th*  information  transfer  on  a 
noisy  channel.  In  the  present  case  the  proouct  aperture  FOV  replaces  the  bandwidth 
time  product  of  Shannon's  formula. 
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Appendix  A 


The  Equation  for  the  Object  Function  Estimator 


H(u,x)  = Hq(u. x)  + tj  Hj(u.x),  ( 

where  Ho(u.x)  is  the  kernel  satisfying  (18-19).  Hj(u.x)  is  an  arbitrary  function, 
and  tj  is  a small  quantity.  By  recalling  (19)  we  obtain: 


because  of  (Al)  and  the  meaning  of  Hq(u,x).  Thus  from  (A2): 

E | ^S(u)  gQ(u)  — S(u)  //  H^u,^)  o(x>  d2x  j //  Hj  (u.{)  o({)  d2{  | . 

Because  Hj(u.x)  is  arbitrary,  (A3)  is  equivalent  to 

* * 

S(u>  E [ go(u>  o (*)j  = S(u)  //  Hq(u.x)  E [o(x)  o ({)  j d2x  . 


By  taking  into  account  (16-17)  we  obtain: 

+ oc  + 0C 

a<?)  — 2 jf  s(u)  e rgo(u)g;  (v)i 

4tt-oo-oc  ^ -J 

-f  oc  + or 

x / / A(x)  HQ(u,x) 


e3-**  d2v=S(u)A(?) 


ir„(x  - *)  (P  + P - P2)  + fi(x  - ? ) <P2  f n* 


,2 

d x , 


(A5) 


or  because  of  (8): 


S(u)  A(£)  Ps  e3  — — 


= S(u)  A(£) 


PS+P1 


If  H0^— 1 1 - + ^Pi  + nn^  H0— — ^ 


Dropping  the  subscript  0.  not  having  any  particular  meaning,  and  interchanging  the 
roles  of  x and J.  Eq.  (20)  is  obtained. 
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Appendix  B 

Some  Properties  of  the  Functions  ^ (u) 

By  using  the  definition  (25)  and  Fourier  Transforming  (21)  we  find: 

X.  *.(u)  = -Kr  ff  rr.(u  - v)  *.(v)  d2v  , (Bl) 

1 1 " 4-  S A " “ 1 “ _ 

an  integral  equation  with  the  same  form  as  (21),  and  having  identical  eigenvalues. 
Thus  the  set  J'J'j(u)  } is  complete  and  orthogonal  in  S.  In  order  to  find  the  normal- 
ization of  <>j(u)  in  S,  consistent  with  (25),  let  us  consider  the  integrated  product: 

//  *.<u)  *k  (u)  d2u  = //  d2u  //  ^.(x)  d2x  ff  ^(|)  e1-- d2f  . (B2) 

S S A A 


By  interchanging  the  order  of  integrations  in  the  left  side,  one  obtains: 


ff  Vj(x)  d2X  ff  rrs(x  -V  vk(i>  • 
A A 


(B3) 


By  invoking  (21-22)  one  concludes  that: 
ff  d2u  = 6jk  • 


(B4) 
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Thus  (25)  and(B4)  show  that  the  functions  ♦ .(u)  form  a double  orthogonal  system, 
a generalization  of  a known  property  of  the  prelate  spheroidal  functions.  Also 
from  (25)  and(B4)  by  using  Parseval's  theorem: 


4 

*.<u) 

A2 

d u 

+ 00 
I 

+ 00 
/ 

|*.<u) 

2 H2 

d u 

- or  - oo 


(B5) 


that  is  the  characteristic  value  A.  is  the  fraction  of  the  total  "power"  associated  with 
* (u).  which  belongs  to  S.  Since  ine  (x)  are  real  it  follows  from  (25)  that: 

3* **c 

*.<-u)  = *.  (u)  . (B6) 

i — i — 

The  eigenvalues  Aj  of  (21),  or  (Bl),  are  positive  numbers  lesser  than  unity.  It  can 
be  shown  that  as  a consequence  of  the  representation  (55)  for  the  aperture  function 
the  following  property  holds: 

OC 

E A.  = . (B7) 

i=G  1 4w 

that  is  the  sum  of  the  eigenvalues— trace  of  the  integral  operator  in  (Bl)-is  conver- 
gent and  is  proportional  to  the  aperture  - FOV  product.  A fundamental  feature  of 
the  numbers  A.  is  their  cut-off  property.  If  they  are  ordered  with  increasing  values 
of  the  index,  their  values  are  close  to  unity  for 

AS 

A.  r 1,  'hen  i < integer  part  of  — ^ » (B8) 

1 4* 


and  close  to  zero  for  greater  i's: 


A. 


l 


0 for  i > integer  part 


(B9> 


This  fundamental  property  has  been  pointed  out  by  Slepian  and  Pollack  for  the  one- 
dimensional case  (that  is,  for  the  prolate  spheroidal  functions)  and  has  been  shown 

jjMfr 

to  be  valid  for  the  general  case  by  Landau. 

*Slepian,  D.  and  Poliak,  H.  P.  (1961)  Prolate  spheroidal  wave  functions,  Fourier 
analysis  and  uncertainty  — 1,  Bell  System  Tech.  J. , 40:43-63. 

**Landau,  II.  j.  (1967)  Necessary  conditions  for  sampling  and  interpolation  of  certain 
entire  functions  Acta  Mathematica,  pp  37-52, 
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Appendix  C 

Power  Associated  with  Each  Observsbla 


The  terms  (34)  are  the  powers  associated  with  the  observables  0{'s.  The  total 
power  received  by  the  aperture  from  the  object  function  can  be  found  to  be  by.  using 
(B7): 

AS 

Total  received  object  power  = Pc  — * , (Cl) 

b 4e 


and  a similar  expression,  with  Pj  replacing  Pg,  holds  for  the  total  received  power 
associated  with  Nj(u).  On  the  other  hand  the  total  angular  noise  power  received 
from  outside  the  FOV  is: 


Total  power  from  outside  FOV  = — - 

4w 


(C2) 


having  here  taken  into  proper  account  the  extent  of  the  visible  space. 

Assuming  correctly  the  noise  Nj  limited  to  visible  space  G,  rather  than  extend- 
ing to  the  entire  wavenumber  plane,  amounts  to  replacing  in  (17)  and  (A5)  the 
symbolic  function  8(x  -^)  with 


= //_ 


jn(x  - «)  j2_.  4*  J1  T *•£ 

a d u = -t—  — i1' 1 ™rr  - 


2* 


X 


FT 
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<tg(£-!> 


(C3) 


We  now  resort  to  qualitative  reasonings.  If  we  consider  a number  of  terms  close 
to  N given  by  (31),  the  functions  ^.(x)  are  slowly  varying  on  the  aperture  (in  the  one 
dimensional  case  they  have  a number  of  zeros  equal  to  their  orders).  On  the  other 
hand  the  function  (C3)  has  a peak  much  sharper  than  that  of  (Tg(x),  given  by  (14), 
under  the  assumption  of  S<<G.  Therefore  no  substantial  error  is  committed  by 
rerlacing  (C3)  by  a delta  function  in  integrals  involving  the  product  of  (C3)  with 
^.(x),  if  i is  not  substantially  greater  than  N.  In  fact,  more  rigorously  it  is  possi- 
ble to  show  that,  for  a given  product  AS: 


lim  ff  ip.(x)  — rrG(x  - J)  - 4w2  6(x  - J) 

s-o  A 1 4r 


0 . 


The  measurement  noise  power  for  each  degree  of  freedom  is  found  to  be 


ff  n(x)  d2x  ff  *M£)  n (£)  d2£ 

J A l“  " A l~ 


= o 


2 

because  of  (16)  and  (22).  Thus  is  equal  to  the  rms  measurement  noise  associ- 
ated with  each  observable  O.,  furnishing  an  alternative  interpretation  of  its  physical 
meaning.  A strictly  related  discussion  by  Yaglom  can  be  found  in  his  book  on  the 
theory  of  stationary  random  functions. 


*Yaglom.  A.  M.  (1965)  An  Introduction  to  the  Theory  of  Stationary  Random  Func- 
tions,  Prentice-Hall,  1*1.  J. , pp  207-213. 


